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1 Introduction 

Based solely on symmetry considerations, Landau theory can provide a reliable de- 
scription of a system's equilibrium behavior near a phase transition. This phenemeno- 
logical approach relates measurable quantities to one another using a minimum set of 
input parameters that can be determined either by comparison to experiment or from 
first-principles approaches. Landau theory (LT) can therefore serve as a conceptual 
bridge between microscopic models and observed macroscopic phenomena. Because 
it assumes spatial averaging of all local fluctuations, LT is particularly well-suited to 
systems with long-range interactions such as superconductors and ferroelectrics. In 
this chapter we hope to give the curious reader a self-contained contemporary pre- 
sentation of the application of Landau theory to ferrolectrics, including references to 
more extensive treatments for those who wish more detail 12 El IH 13 El • We begin 
by developing LT for homogenous bulk ferroelectrics, and then consider the finite- 
size (thin film) case within the related Landau-Ginzburg (LG) approach. Next 
we discuss the treatment of inhomonogeneity within this framework. We end with a 
number of open questions for future pursuits. 

Let us begin by stating in general terms what Landau theory is and then subse- 
quently what it is not. In a nutshell. Landau theory is a symmetry-based analysis 
of equilibrium behavior near a phase transition. In his classic 1937 papers. Lan- 
dau [H] notes that a system cannot change smoothly between two phases of different 
symmetry; the continuous path that a system can take between liquid and gaseous 
states around a first-order critical point is only possible because there is no symmetry 
change involved. Furthermore, because the thermodynamic states of two phases that 
are symmetrically distinct must be the same at their shared transition line, the sym- 
metry of one phase must be higher than that of the other. Landau then characterizes 
the transition in terms of an order parameter, a physical entity that is zero in the 
high-symmetry (disordered) phase, and changes continuously to a finite value once 
the symmetry is lowered; for the case of the paraelectric- ferroelectric transition |^llUj. 



1 




Figure 1: Schematic of relationship of Landau theory and first principles calculations, 
microscopic models and experiment in the vicinity of a phase transition 

this order parameter is the polarization P. The free energy, in the vicinity of the 
transition is then expanded as a power series of the order parameter (JF(P)) where 
only symmetry-compatible terms are retained. The state of the system is then found 
by minimizing this free energy J^{P) with respect to P to obtain Pqi and specific ther- 
modynamic functions are subsequently computed by differentiating JF accordingly. As 
an aside, we note that if the order parameter is multicomponent, then the Landau free 
energy is constructed from all scalar terms consistent with the system's symmetries 
that are powers and scalar products of the order parameter components. Here, for 
simplicity, we'll consider a scalar order parameter (P) unless otherwise specified. 

The appeal of the Landau approach is that it is a straightforward phenemenology 
for linking measurable thermodynamic quantities in the vicinity of a phase transition. 
However it is only as good as its input parameters, the coefficients of the series expan- 
sion terms in J-'{P), which can be determined from experiment or from first principles 
calculations. We emphasize that Landau theory is strictly a macroscopic approach 
and thus it cannot describe any microscopic physics (e.g. atomic displacements etc.) 
associated with the phase transition. In its simplest rendition, the order parameter 
(e.g. P) is considered to be completely uniform. 

A central ansatz of the Landau approach is that the free energy can be represented 
as a series expansion of the order parameter in the vicinity of the transition. At first 
glance, it may be surprising that singular behavior associated with a transition can 
emerge from such a regular expansion. This occurs because the value of the order 
parameter that minimizes the free energy (e.g. Pq) is itself a singular function of the 
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Phenomenology 


Ferroelectric 
(near Tc) 


Landau - Devonshire Theory 
(Uniform Polarization) 


Poled Bulk System 


Landau-Ginzburg Theory 
(Polarization with Spatial Gradient) 


Bulk System 


Landau-Ginzburg Theory 
with Boundary Conditions 


Film 



Table 1: Phenomenologies and Relevant Ferroelectric Systems (close to Tc) 

expansion coefficients which are temperature-dependent. However, more generally we 
expect that this power-law form of T may not be valid very close to the transition. 
Furthermore Landau theory is based on the premise that local fluctuations in the 
order parameter are small, whereas we expect that this will not be the case in the 
immediate vicinity of the transition. 

So when does Landau theory break down? Levanyuk and Ginzburg have devel- 
oped ^21 ^1 a criterion which uses Landau theory to estimate its own demise. More 
specifically, they argue that for LT to be valid, fluctuations of the order parameter 
must be small compared to the order parameter itself when both are coarse-grained 
over a volume determined by the correlation length.. We shall present a detailed ver- 
sion of their argument once we introduce a specific Landau free energy. Qualitatively 
their criterion suggests that LT works well when the coordination number is high, 
which can be achieved either by large dimensionality or by long-range interactions. 
Since significant dipolar interactions are present in ferroelectrics, does this mean that 
Landau theory is always valid in these systems? It is certainly very useful in ferro- 
electrics but we emphasize that the angular dependence of the dipolar forces leads to 
non-Landau behavior very close to the transition, though this region is usually practi- 
cally unobservable (we will be more specific later). Indeed historically the theoretical 
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study of a simplified version of this problem, the uniaxial dipolar model, yielded the 
first predictions of logarithmic corrections to Landau exponents p!F . Detailed mea- 
surements confirmed these quantitative results in appropriate experimental systems 
[TH] . and this set of findings served as an important benchmark in the development 
of the renormalization group method. 

Many readers may have encountered the phenomenology of Landau theory in the 
context of magnetism, and thus it is worth contrasting its application there and in 
ferroelectrics (and we encourage the interested reader to pursue more on this subject 
in the chapter by N. Spaldin). First and perhaps foremost, the couplings in a ferro- 
magnet are predominantly short-range, since interactions between magnetic dipoles 
are essentially negligible. More specifically we note that the interaction energy be- 
tween two magnetic dipoles, each of magnitude one Bohr magneton, at a distance r 
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where aB= 0.05 nm and a ~ ^§75 working in atomic units. By contrast, 

for an electrical dipole of magnitude p = eAas, the dipolar interaction force is 

n2 
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where the parameter A 1 is determined by examining atomic displacements and 
transverse effective charges The ratio of the ferroelectric to the ferromagnetic 
dipolar forces is then of order ^ ^ (137)^, so clearly long-range interactions are more 
significant in ferroelectrics than in generic magnetic systems. Thus LT is more reli- 
able for a broader temperature region in these electric dipolar systems than in their 
magnetic counterparts, though it can provide guiding information in both cases. In 
conventional ferromagnets magnetic dipolar interactions play a role in the determina- 
tion of domain structures, whereas in ferroelectrics the analogous dipolar couplings 
are important at essentially all length-scales. 

Screening represents a second qualitative distinction between ferroelectric and 
ferromagnetic systems. This phenomenon is due to the presence of free charges that 
effectively reduce local fields; it does not have a magnetic analogue due to the ab- 
sence of free magnetic monopoles. Indeed this simple fact about screening has been 
used to explain why magnetism was known in ancient times but ferroelectricity is 
a relatively modern discovery The surface of a ferroelectric crystal attracts free 
charge from its environment to neutralize the effective surface charge that results 
from the polarization discontinuity; this partially accounts for the historical name 
"pyroelectric" to describe transient fields that were observed upon changes in temper- 
ature. Screening effects are particularly pronounced in ultrathin ferroelectric films, 
|17| ITHl [W[ |2ni 1^ 122] where internal electric fields generated by the ferroelectric 
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dipoles are large enough to produce both electronic and ionic currents. These pro- 
cesses contribute to fatigue, to "imprint" and to the memory of domain patterns 
stored over time [Zj. A detailed understanding of these screening processes, partic- 
ularly near metal-ferroelectric interfaces in thin-film devices, requires knowledge of 
metal-induced gap states and charge compensation processes. Since these are local 
effects, they can be captured within a Landau-Ginzburg framework and this physics 
will be discussed later in the chapter. 

The third major difference between LT treatments of ferroelectricity and ferro- 
magnetism arises from the strong coupling between the polarization and the lattice, 
a feature which is usually not very important in the magnetic case. Ferroelectric 
dipoles are generated by the displacement of atoms, and thus on a local scale there 
is strong spatial anistropy. Many magnets retain continuous spin symmetry (e.g. xy 
or Heisenberg), and thus have low barriers associated with the reorientation of the 
magnetization. By contrast in ferroelectrics domain walls, discontinuities in the po- 
larization direction, are roughly a few lattice constants in contrast to the mesoscopic 
length-scale associated with their magnetic counterparts. In general the development 
of a macroscopic polarization in a ferroelectric will be accompanied by a macroscopic 
strain, and thus ferroelectricity and ferroelasticity are closely related phenomena. In- 
deed if this strain is suppressed, the transition temperature and even its order can 
be modified. Elastic degrees of freedom thus must be represented in a LT for most 
ferroelectrics. 

Let us now turn to a general outline of this chapter. We will begin with a discus- 
sion of the homogenous Landau theory for bulk ferroelectrics with spatially-uniform 
polarizations, reviewing first- and second-order transitions and the dielectric response. 
Next we allow the polarization to vary with position within a bulk sample, and dis- 
cuss Landau-Ginzburg theory and correlation functions. Surface boundary conditions, 
corresponding to the situation for thin films, are then added and the resulting polar- 
ization profile is discussed. Here screening effects, depolarization fields and strain will 
be important. The enigmatic extrapolation length, borrowed from the literature on 
surface superconductivity [221, ^e described with remarks about its underlying 
physical meaning. Then we'll switch to a treatment of inhomogenous systems in the 
bulk, where strain effects and domain formation become very important. We end with 
a number of open questions for further exploration. For example, is LT appropriate 
when both long-range strain and Coulomb interactions are present? How can we 
begin to treat the ever-important question of dynamics of domain motion within this 
framework? Landau-Ginzburg theory is certainly only useful when its constructing 
assumptions and its input parameters are appropriate; if so, it is a fairly powerful 
approach, based almost solely on symmetry, to yield important information about 
measurable quantities in the vicinity of a phase transition. In this chapter we em- 
phasize the underlying philosophy and methodology of this phenomenology, pointing 
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the interested reader to detailed discussions of successful applications to particular 
materials including SrTiO^ BaTiO^ and PbZri_xTix03 12^]; towards this 
more pragmatic end, an Appendix of Landau parameter sets for specific common 
ferroelectric systems is included for easy reference. 

2 Landau-Devonshire Theory 
2.1 General Phenomenology 

To our knowledge, Landau's symmetry-based treatment of phase transitions pi IE] was 
first applied to the case of ferroelectrics by Devonshire ^0], and in this section we 
will discuss this Landau-Devonshire theory that is appropriate for bulk systems with 
spatially uniform polarization. In general, the thermodynamic state of any system in 
equilibrium can be completely specified by the values of specific variables; for bulk 
ferroelectrics these include the temperature (T), the polarization (P), the electric field 
(E), the strain (rj), and the stress (a). Usually we apply electric fields E and elastic 
stresses a externally, so we can regard the polarization and the strain as "internal" 
or dependent variables. A fundamental postulate of thermodynamics applied to a 
ferroelectric is that its free energy J-' can be generally expressed as a function of ten 
variables (three components of polarization, six components of the stress tensor, and 
finally one of temperature). Our goal here is to write down an ansatz for this free 
energy, using symmetry arguments whenever possible to pare down the number of 
terms. The second important thermodynamic principle that we will also employ is 
that the values of the dependent variables in thermal equilibrium are obtained at the 
free energy minimum when the free energy is optimized. 

We make the key approximation that in the vicinity of a phase transition we can 
expand the free energy in powers of the dependent variables with coefficients that 
can be fit to experiment or gleaned from microscopic calculations. In the best case, 
we may be able to truncate this series with only a few terms. In order to be more 
specific, let us take a simple example where we expand the free energy in terms of a 
single component of the polarization; for the moment we ignore the strain field, an 
assumption which might be appropriate for a uniaxial ferroelectric. We shall choose 
the origin of energy for the free unpolarized, unstrained crystal to be zero, and hence 
we write 

J^p = lap2 + l5p4 ^ 1 p6 _ 

2 4 

where we have truncated the power series at the sixth term, and a, b and c are coef- 
ficients that will be discussed in more detail shortly. Here we will use the notational 
convention that JF and F = J dVJ-' refer to the free energy density and the total 



6 



free energy respectively. The equilibrium configuration is determined by finding the 
minima of jFp, where we shall have 



. (4, 

This equation gives us an expression for the electric field i^^ as a function of the 
polarization 

E = aP + bP^ + cP^ (5) 

Thus we can determine the linear dielectric susceptibility above the transition by 
differentiating this equation with respect to P and then setting P = to obtain 

In the Landau- Devonshire theory it is assumed that around the Curie point (T ~ Tq) 

a = ao{T-To) (7) 

and the other coefficients in the free energy expansion are independent of temperature. 
Combining the last two equations, we find an expression for the dielectric stiffness 

K = - = ao(T - To) (8) 
X 

which captures the Curie- Weiss behavior (in x) observed in most ferroelectrics for 
T > Tq; this provides additional support for the linear temperature ansatz for a. 
For the sake of completeness, we note that the temperature Tq where a changes sign 
is close but not exactly coincident with the Curie temperature, as we shall discuss 
shortly. 

If we include the linear temperature-dependence of a, we have the general expres- 
sion for the free energy 

J^P = laoiT - To)p2 + hp^ + lcp6 _ EP (9) 
2 4 D 

where, to our knowledge, Oq and c are both positive in all known ferroelectrics Oil]. 
Figure 2 shows the free energy as a function of polarization in the paraelectric (T >> 
To) and the ferroelectric (T << Tq) phases. How this free energy transforms between 
these two configurations will be determined by the sign of the coefficient b; as we shall 
see shortly, its sign will determine the nature of the parelectric-ferroelectric transition, 
and whether the polarization at T < Tq develops continuously or discontinuously. 
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(b) 
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Figure 2: Free energy as a function of polarization for (a) a paraelectric material, 
and for (b) a ferroelectric material 

2.2 Second-Order (Continuous) Transition 

If 6 > 0, then a second-order transition occurs at T = Tq, and the free energy will 
evolve continuously as a function of decreasing temperature from the first schematic 
(P = 0) in Fig 2 to the second, that has minima at finite polarizations P = ±Po- 
The spontaneous polarization can be estimated by setting = in (0); since all the 
coefficients are positive, we will only retain the two lowest-order terms. The result is 



where we see that the spontataneous polarization Pg will increase with decreasing 
temperature from the point T = Tq. We note that if we determine the dielectric 
stiffness below the transition (T < Tq) then we find 



which is to be compared with (jSJ, its value just above Tq; perusal of these two 
expressions suggests that k vanishes at T = Tq and that consequently the dielectric 
suseptibility diverges. In practice, the latter achieves large but finite values in real 
materials where this expression is appropriate j3| H] . We can now also solve for the 
discontinuity in the specific heat at the transition; using P = for T > Tq and (fTUI) 
for T < To, we determine 
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K = 2ao(To - T) 



(11) 



AC, = C,{T = T+) - C,{T = T,) 



(12) 



where C, 



v — 



T to obtain 



AC, 



(13) 
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Figure 3: Second order phase transition, (a) Free energy as a function of the polar- 
ization at T > To, T = To, and T < To; (b) Spontaneous polarization PoiT) as a 
function of temperature (c) The susceptibility x and its inverse, where x = 9P/dE\p^ 
is evaluated at the equilibrium polarization Po{T) 

for the Landau expression for this quantity. Schematics associated with this second- 
order phase transition as described by Landau-Devonshire theory are displayed in 
Fig. El 

2.3 First-Order (Discontinuous) Transition 

Logically (and practically as it turns out), we should also consider the case of 6 < 
(while c remains positive). This is sketched in Fig. |3] 

With the negative quartic coefficient, it should be clear that even if T > Tq (such 
that the quadratic coefficent is positive) the free energy may have a subsidiary mini- 
mum at non-zero P. As a is reduced, which corresponds physically to the reduction 
of the temperature, this minimum will drop in energy below that of the unpolarized 
state, and so will be the thermodynamically favored configuration. The temperature 
at which this happens is, by definition, the Curie temperature T^ which however now 
exceeds To. At any temperature between Tc and Tq the unpolarized phase exists as 
a local minimum of the free energy. The most important feature of this phase tran- 
sition is that the order parameter jumps discontinuously to zero at T^. This type of 
phase transition is usually called a first-order or discontinuous transition, of which 
solid-liquid transitions are common examples. 

The procedure for finding the spontaneous polarization and the linear dielectric 
susceptibility is conceptually the same as before, but now one cannot neglect the 
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T>Tc . (a) 




Figure 4: First order phase transition, a) Free energy as a function of the polarization 
at T > Tc, T = Tc, and T = Tq < T^, (b) Spontaneous polarization Po(T) as a function 
of temperature (c) Susceptibility x- 

quartic and sixth order terms. Qualitatively we then find a dielectric stiffness (inverse 
of the linear susceptibility) that does not vanish at Tq, corresponding to the finite 
jump in both the susceptibility and the spontaneous polarization at the transition. 
Schematics of the free energy, the spontaneous polarization, the dielectric stiffness 
and the linear susceptibility are shown in Fig. |3] We note that at T = T,, the three 
minima are energetically degenerate. As a result, the system's behavior at T = Tc 
will depend on whether it is approaching Tc from lower or higher temperatures. More 
specifically, the system will be in one of the two finite polarization (P 7^ 0) minima if 
it is heated from an initial low temperature Tj < Tc, whereas it will be in a paraelectric 
state (P = 0) if the initial temperature is high (Tj > Tc). Indeed the phenomenon of 
thermal hysteresis, where the transition temperature depends on whether the sample 
is heated or cooled, is prevalent in a number of first-order ferroelectrics including 
barium titanate |3,. We emphasize that it is only for Tq < Tc that the ferroelectric 
minima are thermodynamically favorable. 

In a ferroelectric below Tq there are (at least) two minima of the free energy, cor- 
responding to spontaneous polarizations of different spatial orientations. The barrier 
between these minima means that a small electric field will not immediately switch 
the polarization. We note that the Landau-Devonshire theory described here predicts 
hysteresis jSl EI, shown schematically in Figure 5, in the ideal (and fictitious) case 
where all the dipoles have to be overturned together to switch from one polarization 
orientation to the other. 

We note that for ferroelectrics with either continuous or discontinuous transitions, 
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Figure 5: Schematic picture of hysteresis in an ideahzed ferroelectric 

values of the Landau coefficients can be determined by comparing experimental mea- 
surements of K and Ps to the Landau-Devonshire expressions; such parameter sets for 
a number of commonly studied materials are included in an Appendix at the end of 
this chapter. 

2.4 Coupling to strain 

An important feature of ferroelectric materials is their great sensitivity to elastic 
stress. In order to understand why this is so, we can again take recourse to Landau 
theory by adding in strain dependent terms to ©• The strain in a solid is measured 
by how the displacement m of a point in the solid varies with position r, and since 
this is the dependence of a vector upon a vector, the answer is a tensor: the strain is 
usually defined as 



where here i,j indicate the x,y,z components of the vectors, rj is therefore a 3x3 
symmetric matrix, with six independent components. In materials that are cubic 
(or nearly so) there will be three independent components to the strain: the volume 
strain (uniform in all three directions and two kinds of shear). 

In general, the polarization will couple to one or more types of strain, and specif- 




(14) 
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ically which ones can generally be seen by inspection. Consider a cubic crystal (e.g. 
BaTiO^) that undergoes a ferroelectric phase transition to a state where the polar- 
ization can point along one of the six orthogonal cubic directions. Now it is clear that 
there is a special axis (one of the six directions after the symmetry has been broken) 
and so it would no longer be expected that the crystal as a whole will remain cubic 
— one expects a distortion into a tetragonal crystal, which can be described by a 
tetragonal strain r]. The fact that the lowest order coupling allowed in this case is of 
the form rjP^ (and not, for example, rjP or rj'^P) is of course a matter of symmetry. 
For a uniaxial ferroelectric, the leading order terms will be of the following form 

J^^ = ^Kr]^ + Qr^P^ + ..-r]a , (15) 

Here t] is (a component of) the strain field, and the first term represents Hooke's law, 
namely that the elastic energy stored in a solid is quadratically dependent on the 
distortion, so K is (one of) the elastic constant (s). The second term is a coupling 
between the elastic strain and the polarization; the fact that this is linear in the strain 
and quadratic in the polarization depends on the special symmetry of the transition. 
While this is the leading term in pseudocubic materials, there are other materials 
(e.g. KH2PO4) where the symmetry is lower, and the coupling can be of the form 
r)P - linear in both strain and polarization. Materials with a linear relation between 
stress and polarization are called piezoelectric. 

Using the free energy, which now consists of the terms in Q and (fT3j) . JF = jFp+jF,. 
we can now determine the properties in equilibrium by minimizing with respect to 
both P and 77, viz 

dHP, V) d^P, V) . .... 
^i^ = ^^ = ° ^^'^ 

Let us take the second of these equations first: 

?^f^ = K„ + QP^-a (17) 
or) 

There are a few different limits to look at. Firstly, note that if the polarization is 
zero, we get Hooke's law t] = a/K. The second - apparently trivial case - is when 
a stress is applied to force the strain to be exactly zero at all times. This is not as 
absurd as it seems, because often crystals can be considered to be clamped by their 
surroundings so that no strain is allowed at all. One common situation is of a thin 
epitaxial film which is forced to have the lattice constants matched to the substrate, 
and is free to relax only in the third direction. In the case of perfect "clamping" 
?7 = 0, and the free energy is just as before. 
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The third case to consider is when no external stresses are apphed (a = 0), and 
we then have 



so that a spontaneous (tetragonal strain) occurs proportional to the square of the 
polarization. Notice now that we can substitute for the strain as a function of polar- 
ization, and we have a free energy 

T{P,r^{P)) = ]-aP^ + -{h-2Q'/K)P' + \cP' + ...~EP . (19) 
2 4 o 

In comparison with the clamped system, the only change is to reduce the quartic 
coefficient (notice that the result is independent of the sign of Q). This means that 
in the case of an already first-order transition (6 < 0) the transition is driven even 
more strongly first order, and Tq is raised. In many ferroelectrics, the effects of 
clamping can be large - shifting the transition temperature by tens of Kelvin, and 
even changing the order of the transition. If 2Q'^ / K > 6 > 0, a first-order transition 
becomes second-order in a clamped system where the strain is allowed to relax; this 
is the case jS] in BaTiO^. 

Application of external hydrostatic pressure leads to a shift in Tq whose sign is 
determined by that of the volume expansion coefficient; a reduced lattice constant is 
accompanied by a reduction in Tq, consistent with measurements where similar effects 
have been achieved by chemical substitution j27] . The effect of biaxial stress, where 
the crystal contracts along two axes and expands along the third, on Tq was also 
studied both theoretically and experimentally j2H]; enhancement of Tq was predicted 
and observed in barium titanate crystals [SHj. We shall see shortly that this effect is 
even more dramatic in thin films where such elastic effects occur at the interface to 
the substrate. 



2.5 Domains 

So far we have pretended that the polarization in a ferroelectric can be treated as 
entirely uniform, and this is far from the case. There are many reasons for the exis- 
tence of domains, including non-uniform strain, microscopic defects, and the thermal 
and electrical history of the sample. But even in an ideal crystal, domains are to be 
expected for energetic reasons associated with electrostatics. 

The macroscopic bulk polarization is produced by a displacement of positive 
charge relative to negative charge; at the surface of the sample there must then 
be a net charge density of opposite signs on the opposite sides of the crystal (Fig. E)). 
Any inhomogeneity in the polarization acts like a charge density - and in particular 
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Figure 6: Surface charge density generated by a bulk polarization at an interface 



the discontinuity in the polarization at an interface or boundary acts like a surface 
charge, as can be seen by the following elementary electrostatic argument. 
The potential V{r) induced by a dipole p at the origin is 



V{f) 



Arcen \r 



(20) 



and thus the potential due to a distributed polarization P(r) in a volume r bounded 
by a surface S is 



V{r] 



^ ^ dr' Pi?) ■ Vr' I — - 



Using the product rule: V ■ (A/) = /V ■ A \ A - V/, we have 



(21) 



V{f) 



1 



1 



dr' 



V • Pfr' 



V 



^^V-P(r') 
dr'^ — ^ + 



dS ■ P{? 
s \r — r'l 



(22) 
(23) 



where the last manipulation uses the divergence theorem leading to a surface integral. 
The potential appears to be generated by a distributed bulk charge density 



p(f) 



-V ■ P(r) 



as well surface charge density 



a = P ■ n 



(24) 



(25) 



where h is the vector normal to the boundary surface S. Occasionally, some books 
that like to make arbitrary distinctions between "bound" and "free" charge will call 
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the charges generated by an inhomogeneous polarization "fictitious" . These charges 
are fictitious only in the sense that they are conserved - namely if we begin with 
an electrically neutral system and create a polarization density from it by moving 
the electrons relative to the ions, the integral of the polarization charge (over the 
whole system) must vanish. They are in all other respects quite real. The surface 
charges generate electric fields, both internal and external to the sample, and the 
fields themselves store energy - just as in a capacitor. In an approximation where 

— * 

the polarization P is assumed to be homogeneous within a domain, and changing 
abruptly at surfaces or interfaces, only the surface term survives. 

The system will minimize its energy by eliminating - as far as is possible - the 
surface charges, and in a thin film for example this makes clear that the preferred 
orientation of the polarization will be in the plane of the film, rather than pointing 
perpendicular to the film. If one has a crystal that is thick in all dimensions, another 
situation is preferable - to introduce domains where the polarization is oriented to be 
always parallel to the crystal surface. Such a situation can of course only be achieved 
by introducing domain walls into the bulk, and this will produce polarization charges 
unless the walls are appropriately oriented. 

The interface charge density between two neighbouring domains is, by extension 
of the formula at a free surface, 

a = (A - P2) ■ ft (26) 

where Pi, P2 are the polarizations of the two domains, and ri is a unit vector normal 
to the interface. The two cases when the surface charge will exactly vanish are when 
the polarizations are antiparallel to each other, and parallel to the domain wall (called 
a 180° domain wall), or when the domain wall bisects the angle between two domains 
pointing head-to-tail. Because we are often dealing with nearly cubic crystals where 
the possible polarizations are at 90° to each other, the latter is often termed a 90° 
domain wall. Both of these wall types are illustrated in the sketch Fig. [71 which is 
the ideal configuration of the polarization in a " bar" of ferroelectric single crystal — 
and which bears obvious comparison to the ideal magnetic configuration in a single 
crystal bar magnet, that may be familiar to some. 

The presence of domain walls involves other energy costs. Firstly, the domain 
wall is microscopically different from the bulk, and the energy gain of forming the 
polarized ground state has been lost. Secondly, the polarization is coupled to elastic 
strain - as we saw above - and we must also make sure that the strain fields are 
compatible. There is no difficulty with a 180° domain wall, but a 90° domain wall 
produces problems, and the "ideal" configuration of Fig. [7| is by no means ideal 
in terms of the strain fields. The existence of domains in a sample automatically 
generates inhomogeneities in the strain which interact in complex ways. 
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Figure 7: Ideal domain configuration in a single crystal of cubic ferroelectric material, 
where the coupling to strain is negligible. On the right is the configuration adopted 
when strain effects are important. 

Some contemplation of the two-dimensional pictures above makes it evident that 
only certain kinds of domain structures are allowed, unless dislocations, cracks, or 
voids are introduced into the crystal. The moment we have inhomogeneous strain, 
we must remember to maintain the condition that the displacement field u{r) must 
be continuous; however, the free energy depends only on the strains rj, which are 
derivatives of u (see Eq. ()14|1 ). Imposing the continuity constraint gives rise to the 
additional (St Venant) compatibility conditions on the strain fields (23 

V X (V X r^) = 0, (27) 

which formally reduces the number of independent strain components. When dealing 
with the allowed arrangements of domains of (locally) homogeneous strain (except 
at the boundary), the application of condition Eq. (j27j) (typically by taking a line- 
integral around the intersection of three domains) determines the legality of various 
domain arrangements in two and three dimensions. 

This condition can also be used to integrate out a single component of the strain, 
at the expense of introducing long-range forces jHOI, and this approach can then be 
used to support calculations and simulations of domain structures in two |31jJ and 
three ^32j dimensions. This methodology is now quite far advanced in understanding 
ferroelastic effects on domain structures in martensites jHSl though it has not yet been 
combined with simulations of both the elastic and the ferroelectric order parameter. 

There is another feature of ferroelectric domain walls which is not shared by 
magnetic domains. As far as we are aware, magnetic monopoles do not exist, and 
therefore the (fictitious in the sense of conserved) magnetic charges that are generated 
on magnetic boundaries cannot be screened. In ferroelectrics, the surface polarization 
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charges can be, and are, screened by real electrical charges - from impurities, defects, 
and migrating ions for example. The motion of charged species (on a fairly slow time 
scale) to domain walls then provides a mechanism for memory and associated physical 
phenomena. 

This brief introduction to domain structure highlights several issues. The first is 
that we need to extend the Landau-Devonshire theory of a uniform ferroelectric to 
incorporate non-uniform polarization. This we do by a parameterization of the free 
energy now including terms in powers of gradients - this is Ginzburg-Landau theory, 
introduced in the next section. Secondly, we noted that the geometry of the sample 
influences the domain structure, and in particular that there is potentially antagonism 
between the strain fields and electric fields introduced by the domain structure. The 
final major point is to note that the effect of inhomogeneities in ferroelectrics is to 
introduce long-range forces, from either charge inhomogeneity or strain fields. It 
is these two phenomena that complicate a local description of the ferroelectric free 
energy. 

3 Landau-Ginzburg Theory 
3.1 General Considerations 

The Landau-Devonshire theory that we discussed in the previous Section is well-suited 
for the description of a poled bulk ferroelectric near its transition (To) with a polariza- 
tion for T < Tq that is spatially uniform. However we also offered physical motivation 
for treating spatially nonuniform polarization in unpoled bulk ferroelectrics, and this 
will be the subject of our treatment in this Section. Generally speaking, Landau- 
Ginzburg (LG) theory incorporates small spatial variations of the order parameter 
(here the polarization) within the phenomenological Landau-Devonshire theory with 
assumptions used originally by Ornstein-Zernike; there fluctuations of the order pa- 
rameter at different wavevectors are assumed uncorrelated and thus independent. 

Let us now discuss specifics of this Landau-Ginzburg approach. Slow variations in 
the direction of the polarization lead to an additional contribution in the free energy 
density (AJ^lg) that is proportional to |VPp. Though here, as before, we restrict 
our discussion to the case of a single- component order parameter, we note that this 
specific quadratic form of AJ-'cg is valid for more general vector order parameters 
based on symmetry considerations |2] . Let us return to the scalar order parameter. 
Then if we define a polarization density P{r) where r is a dimensional spatial vector, 
then to lowest order the Landau-Ginzburg free energy is 

F = ao{T - To) J d''f{P{r)f + rf¥{VP(f)}' (28) 
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where the second term is simply the leading contribution in an expansion of polarization- 
polarization interactions that estimates the additional free energy cost if the polariza- 
tions at different spatial positions are not parallel. Using the fact that the polarization 
is real (i.e. P*{q) = P{—q)), then we can perform a standard Fourier transform to 
rewrite the Landau-Ginzburg free energy as 

F = J ^(ao(T-To)+7fc^)|P(fc)r- (29) 

We note that here we have two quadratic degrees of freedom; by the equipartition 
theorem we can then write 

{ao{T-To)+^e)\P{k)\^ = knT (30) 

where fc^ is the Boltzmann constant. 

3.2 The Polarization Correlation Function 

We are now in a position to determine the static two-site polarization correlation 
function 

gir) = (P(r-)P(0)) - (P(0))2 (31) 

which is a measure of the fluctuations of the polarization. We note that in ferro- 
electrics such quantities can be measured with diffuse X-ray scattering jHOlj. Taking 
a Fourier transform of (y'(r) above in the paraelectric phase (where (P(0)) = 0) and 
assuming that the different wavevector modes are uncorrelated, we obtain 

9{<l} = {\Pm') (32) 
which then, combined with the expression in ()30p . yields 

giq) = ^^"^ (33) 

ao{T - To) + 'yk^ 

as the expression for the Fourier transform of the correlation function. Now we simply 
take the inverse transform to find 

9ir) - -^-^ (34) 

and 

- T = To (35) 

fy lyU, Zi 
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where the correlation length is 



'''' 

and t = is the reduced temperature. Physically for T > Tq this correlation 

length, ^, corresponds to the length-scale over which the polarization exists; it diverges 
at T = Tq, and for T < Tq it is associated with the length-scale over which P varies 
from its equilibrium value. We note that the polarization correlation function is 
exponentially decreasing above the transition, but has a dimensionally-dependent 
power-law form when T ~ Tq. As a final note to this subsection, we remark that the 
dimensions of 7 and ao{T — Tq) are -jMj^ and respectively, so that [^0] = = L 
as expected. 



3.3 The Levanyuk-Ginzburg Criterion 

How reliable is this Landau- Ginzburg theory, and when is it no longer valid? The LG 
approach is a long-wavelength description of a system near a phase transition, where 
its modes are coarse-grained on scales of order the correlation length ^. The resulting 
effective free energy density is written as an expansion of the order parameter (e.g. P) 
averaged over a volume determined by ^. When are the fluctuations of the order 
parameter, averaged over the correlation volume small in comparison with its 
coarse-grained average? The answer, relevant for the validity of this long- wavelength 
approach, is dependent on system dimensionality. We have just seen that the LG 
approach yields a two-site correlation function, a measure of the order parameter 
fluctuations, that increases dramatically at T ~ Tq; it is therefore clear that in the 
immediate vicinity of Tq the LG expansion is no longer valid. So how close can we 
get to the transition itself before this approach no longer works? 

As it turns out, we can use Landau- Ginzburg theory itself to determine its own fal- 
libility f2llEI- Before we present this argument, let us return to Landau-Devonshire 
theory and look at its predictions for behavior near the phase transition. We recall 
the expression for the free energy 

:Fp = ^ao{T - To)P^ + hp^ - EP (37) 

which we have truncated at quartic order. The equilibrium value of the polarization 
corresponds to a free energy minimum 

--^ = ao{T-To)P + bP^ = (38) 



19 



which yields (for T < Tq) 



P = f'^\^^ - {-tf". (39) 

whereas P = in the paraelectric phase. The form of the free energy, (jHTj) . combined 
with these results for the polarization imply that the associated specific heat, = 
— T|^, has a discontinuity at the transition so that the mean- field exponent a is zero 

(i.e. a ~ |tn. 

Equipped with this information, we are now ready to determine the condition for 
when Landau-Ginzburg theory breaks down. Basically this occurs when the fluctua- 
tion free energy associated with a typical fluctuation of order the correlation length 
is comparable to the total free energy. The fluctuation free energy per unit volume 
can be estimated as 

kT 

:Ffi^ct ~ ^ ~ (40) 

where we have used ^ ~ \t\~'^- We have defined the specific heat, C„ ~ li^^l"" so that 
two integrations yield ~ Thus for a consistent theory, for t — 0, we must 

have 

d* > ^— ^. (41) 

Inputting the mean- field values (a = and z/ = 1/2), we get 

d* > 4 (42) 

which indicates that above the upper critical dimension, d* = 4, the LG approach is 
valid. We note that the value of d* depends on the form of the Landau expansion 
and can be different for systems with different underlying symmetries j35j . 

This calculation indicates that for real materials in < 4, the Landau-Ginzburg 
approach breaks down close to the phase transition... but how close? Intuitively it 
seems that the larger the number of neighbors the better it does, which suggests that 
it works better for systems with long-range forces. Let us use what we have learned 
to put this speculation on firmer footing, and here we follow the original Levanyuk- 
Ginzburg reasoning [T21 CSI • For T < Tq let us argue that coarse-grained fluctuations 
in the polarization must be small in comparison with the average polarization itself; 
this translates into the condition 

g{r) « P2 (43) 

where |r| ~ ^. Using our previous results in (jMj) with d = 3 and (fH^. we can rewrite 
this expression as 

Ml « - (44) 
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which, when using specific expressions for ^(T) and in and ()13|). results in 
the standard form of the Levanyuk-Ginzburg criterion 



where ^ is a constant that is unimportant for our present purposes; here the key 
point is that the range of the interaction, ^o, plays an important role in determining 
the validity of the mean-field theory. More specifically, the Levanyuk-Ginzburg tem- 
perature, determined by the expression in ()45|1 . is inversely proportional to the sixth 
power of the range of the interaction, and thus Landau- Ginzburg theory should be 
quite reliable in the vicinity of the paraelectric-ferroelectric transition. We note that 
if we generalize the condition to d-dimensions then, using (fH^ with |r| ~ ^ and 
(jnni)) we obtain 

B\t\— « \t\^ (46) 



which is satisfied for d > d* = 4 for arbitrarily small \t\ where i3 is a constant. We 
have thus recovered the same result that we obtained earlier in this section, namely 
that Landau- Ginzburg theory is exact for dimensions greater d > 4 for the ferro- 
electrics under consideration. Here we remark that the angular-dependence of the 
dipolar interactions has not been considered here, and it is indeed this feature that 
leads to logarithmic corrections to mean-field exponents in certain experimentally 
observable (i.e. d = 3) cases jT3]. We note that here we have referred to this condi- 
tion as the Levanyuk-Ginzburg criterion (rather than simply that due to Ginzburg) 
since we have learned [36j that Levanyuk, a student of Ginzburg's, derived this con- 
dition independently of his advisor and indeed published the result in a sole author 
publication [T2] a year before Ginzburg's paper appeared |13] . 



3.4 Displacive and order-disorder transitions 

The Levanyuk-Ginzburg criterion indicates the validity of the mean-field approach to 
the transition, but more fundamental in many practical senses is the validity of the 
continuum approximation underlying the LG theory. A continuum theory is of course 
good close enough to a second order critical point since the diverging correlation 
length ^ is much larger than the interatomic spacing — this is the condition that 
coarse graining works. But far from the transition — either well above Tc or well 
below it — this condition will not necessarily hold. 

If one takes the Landau theory literally as a description of the phase transition and 
its dynamics, the order parameter vanishes uniformly everywhere above the transition 
temperature. Equivalently one would say that the lattice displacements corresponding 
to the electronic polarization are equally uniform - and zero above T^, finite below. 
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Such a type of transition is termed displacive. A contrasting limit for an ordering 
transition is that of local moment magnetism j^, wherein the magnetic moment 
on an atom persists above Tcurie but the long-range order is destroyed. In the case 
of a ferroelectric (or other structural transition), this limit is described as an order- 
disorder transition, where the local symmetry-breaking distortions are present in every 
unit cell above T^, but are randomly oriented at high temperatures, so there is no 
net polarization. Any real material will of course lie somewhere along the spectrum 
between these two limits. 

We have argued above that one reason for the practical usefulness of the LG 
approach is that there are intrinsic interactions (elastic and charge) that are long- 
range, and these induce a long correlation length. But there are also conditions 
related to the underlying microscopic drivers for ferroelectricity that are ineluctably 
electronic in nature. Although in this chapter we are largely avoiding microscopic 
theory here we engage in a short digression on the topic of electronic ferroelectricity, 
which has some relevance to the issue. 

Spontaneous development of a ferroelectric moment is always associated with a 
(broken symmetry) atomic displacement, and often this is pictured as an essentially 
rigid displacement of charge. But of course the driving force for ferroelectricity arises 
from interactions in the electron system, so a completely rigid displacement of ions is 
not to be expected. One may study this even in linear response theory - how much 
charge redistribution is produced by a small displacement of the ion - or equiva- 
lently how large an electrical dipole is produced by a particular phonon displacement. 
Phonons in solids can be separated into acoustic and optical branches, where the 
eigenmodes of the latter describe the relative motion of different ionic species against 
each other. One may quantify the electronic rigidity for an optic mode by a number 
known as the transverse electric charge jHj, which measures (to linear order) 
the average electrical dipole moment per unit cell p generated by the corresponding 
relative ionic displacement u (taken to be the same in every cell), viz. 

p = Z^eu (47) 

where e is the electron charge. There is a different value of effective charge Z^ 
for each optical mode ( and it may be zero if the displacement preserves inversion 
symmetry); even a non-ferro-electric crystal will usually have dipole-active modes. A 
simple case is just that of the rocksalt structure, where the optical eigenmode is a 
displacement of the two sublattices relative to each other. The effective charge can 
be straightforwardly measured by the optical response of an infra-red active phonon. 
The connection to ferroelectricity is that if such a phonon goes "soft" and acquires 
a finite frozen amplitude Uq, then Z^euQ is the approximate magnitude (to linear 
order) of the ferroelectric moment per unit cell. This is the classic soft mode theory 
of ferroelectric instability PH} IH^ liUj. 
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In the context of the Landau theory, this relates the polarization to the atomic 
motion. In a rigid ion picture, one expects to reflect the charge on the ion, but 
it is in fact not uncommon to measure effective charges that are quite large: for 
example IH] in the rocksalt structure IV- VI compounds {Z^ values in brackets) PbS 
(4.8), PbSe (5.8), PbTe (6.5), SnTe(8.1), GeTe(ll), where a rigid ion viewpoint would 
lead one to expect Z^ < 2. Only the last two of these materials (which are the least 
ionic in the Pauling sense!) become ferroelectric at low temperature. What these 
large numbers represent physically is of course that as the nuclei are displaced the 
electron distribution sloshes from one side of the ion to another, with the electrons 
moving much further than the nuclei themselves. This large amplification of motion 
happens in electronically driven transitions of the Peierls type j42j; it can also happen 
in cases of electronically-driven spontaneous symmetry breaking where the Coulomb 
interaction favours the breaking of spontaneous orbital symmetry I44j . 

In these cases, the large induced polarization arises because of a broken orbital 
symmetry. Whatever the underlying model, in what one can see in the not entirely 
formal limit of a ferroelectric that has Z^ oo, the transition is indeed electronic 
in nature. Moreover it corresponds in that case to a metal-insulator transition as 
well, where the distortion leads to the opening of a small electronic energy gap 2A 
in the spectrum. Such transitions are in fact charge density wave instabilities in 
disguise, where the broken symmetry happens also to break inversion symmetry, and 
the CDW onset is also that of ferroelectricity. The driving force for these transitions 
is the opening (or increase of) an electronic gap at the fermi energy (for a review, see 

Suppose for a moment that we can ignore the lattice dynamical effects on this 
transition, we can then estimate the transition temperature and correlation length 
in a weak coupling limit, which is well-known to be mathematically just like the 
Bardeen-Cooper-Schrieffer theory of superconductivity [IHl EZl EH] : 

^ oc A ~ We-^/^ (48) 
a A 

where W is the electronic bandwidth, a the lattice constant, and A ^ 1 the effective 
dimensionless interaction strength (scaled by the electronic bandwidth), which must 
be small for this kind of theory to be appropriate. In weak coupling, the gap and 
transition temperature are small (in comparison to the electronic bandwidth) and the 
correlation length is long. 

Now we consider the effect of the coupling to the lattice, and in particular we must 
take account of dynamical phonons. Provided ksT^^^ is much smaller than a char- 
acteristic phonon energy huD, there will be no appreciable fluctuations of dynamical 
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modes until one is at a temperature very close to the transition, so the mean-field 
theory is rather good as an estimate of the transition temperature. And moreover, 
even above Tc the dynamical modes are weakly excited so there are few thermal fiuc- 
tuations of the lattice. All that remains is just the (small, in this case) mean-field 
lattice distortion that smoothly vanishes at Tc. We see that the condition for the 
displacive limit is that of ultra-weak coupling, namely that 



W 



> e-i/-^ . (50) 



However, if this condition is violated, thermal fiuctuations of dynamical modes will 
become important, and the extra entropy associated with the lattice distortions will 
drive a transition to the paraelectric phase at a temperature well below the mean-field 
one . Here we now crossover to the "order-disorder" regime of the transition, and 
we remark that (as commonly observed in CDW transitions) this crossover happens 
when A is still very small, since typically the electronic bandwidth is two to three 
orders of magnitude larger than a typical phonon energy. 

Thus we see that there are two potentially quite distinct reasons to stabilise the 
displacive picture. One is the coupling to elastic strain, and the presence of long- 
range forces, about which we will have more to say later. This is important in all 
classes of ferroelectrics. A second type is where the transition is predominantly driven 
by electronic redistribution of states near the fermi level of a metal or narrow-gap 
semiconductor, and furthermore the displacive limit obtains only when the coupling 
is very weak, so that electronic gaps are small in comparison to the Debye frequency. 
These CDW-like systems are rare as classic ferroelectrics, (except maybe for SnTe, 
which has a transition near 40K and shows the classic soft-mode behavior |50j), but 
they are probably quite common in small-moment ferroelectrics where the ferroelec- 
tricity is an "accidental" adjunct to a charge-ordering transition. Such is potentially 
the case in charge-ordered manganites, for example [SI]. 



3.5 Recent Developments in Bulk Ferroelectricity 

Landau-Devonshire theory has been successful at reproducing the observed phase 
behavior in simple ferroelectric compounds (e.g. BaTiO^ and PbTiO^) and these 
results are described in many textbooks HI E]- Recently it was shown that 
the standard sixth-order free energy expansion can not account for an observed phase 
where the polarization is not symmetry-restricted to an axis. More specifically, a 
monoclinic phase was observed in PZT for a given temperature and compositional 
range [HHl; that did not emerge from a standard Landau-Devonshire (sixth-order) 
treatment [M]. However this phase did emerge in the observed parameter regime 
from simulations ^55j based on a first-principles effective Hamiltonian approach 156] : 
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this result suggested that no additional physics (or specifically instabilities) were 
necessary for its description. An eighth-order expansion of the free energy yielded 
the observed monoclinic phase where it was argued that the higher-order terms 
might be generated by disorder-averaging. This analysis will be relevant for other 
ferroelectric phases where the polarization is constrained to a symmetry plane. 

In the Introduction we alluded to a link between first-principles calculations and 
Landau theory, and here we discuss it more specifically. In the last two decades, there 
has been great progress in the collective understanding of atomic-scale ferroelectric- 
ity through detailed density-functional theory investigations. We refer the interested 
reader to two recent reviews |K7| EH] and other chapters in this book for more infor- 
mation about these developments. First principles effective Hamiltonians, based on 
these first- principles approaches with reduced degrees of freedom, have been analyzed 
by statistical methods to explore the finite-temperature behavior of ferroelectrics. In 
particular such an effective Hamiltonian approach was very successful in reproducing 
the known phase behavior [321 of BaTiO^. More recently the compatibility of this par- 
ticular Hamiltonian and the Landau-Devonshire theory of BaTiO^ was demonstrated 
[nO] in a careful Monte Carlo study of the order-parameter configuration space. Thus 
first-principles density functional calculations, via an effective Hamiltonian, can pro- 
vide the input coefficients for Landau- Ginzburg theory; this bridging of theoretical 
methods is crucial for studying phenomena on many length-scales in ferroelectrics. 

In this Chapter, we focus primarily on the thermodynamic description of ferro- 
electrics using Landau- Ginzburg theory. However the resulting free energy can be 
incorporated into a treatment of thermally activated dynamical processes, specifi- 
cally the determination of nucleation rates fHI]. We note in passing that the coupling 
between the polarization and the long-range elastic degrees of freedom in bulk mate- 
rials implies that nucleation of new ferroelectric regions must be a cooperative effect 
in that these events cannot occur completely independently. Dynamical studies of 
field- quenched BaTiOs crystals provide experimental support for this statement [U?]. 
This observed behavior can be qualitatively described by a modified time-dependent 
Landau analysis with simple feedback to account for the strain-mediation between 
growing regions [63]. 

In concluding our discussion on bulk ferroelectricity, we note that there was much 
work done in this field in the former Soviet Union that was not communicated to the 
West. This situation was recently summarized by Ginzburg [HI], and we point the 
interested reader to these thought-provoking reminiscences. 
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4 Reduced Size and Other Boundary Effects 



4.1 General Discussion 

Ferroelectric materials are very sensitive to electromechanical boundary conditions 
due to the long-range nature of their underlying electrostatic interactions and to 
the strong coupling between the polarization and the strain EZl EH]- Thus the 
influence of surfaces and other boundary conditions on their collective behavior is 
very pronounced Furthermore there are numerous ways in which the external 
environment, via applied boundary conditions, can alter ferroelectric behavior; ex- 
amples include electric flelds due to surface charges and homogeneous misflt strains. 
Finally there is a signiflcant commercial impetus to use ferroelectrics for portable 
high-density data storage; for such increasingly miniaturized applications, key size- 
dependent effects must be understood and modeled in order to optimize design 
Typically the fabricated structures will be of linear dimension 10 nanometer - 1 mi- 
cron, length-scales that are not accessible to purely atomistic methods. However such 
approaches can be used towards determining the effects of local variations (e.g. com- 
positions, strains, displacements) in conjunction with speciflc boundary conditions; 
the results could then be incorporated into a bridging phenomenological theory that 
spans physics on longer length-scales to enable direct comparison with experimental 
observation. 

In this section we summarize key features of Landau- Ginzburg approaches to 
ferroelectrics with speciflc boundary conditions, focusing on the well-studied case of 
planar geometries. We begin with semi-inflnite ferroelectrics with a free surface, a 
case that can be easily generalized to a free-standing thin fllm. Because of reduced 
coordination number at a free surface compared to the bulk, we expect the average 
polarization to display altered behavior at the boundary which could lead to changes 
in its overall thermodynamic behavior jHH inH|- We observe how this need for a 
speciflc boundary condition emerges technically from a Landau-Ginzburg approach 
of the bulk and the surface free energies. We also discuss the necessary assumptions 
underlying the emergence of the so-called extrapolation length from such a treatment 
and suggest a test for this often-used (but rarely justifled) approach [ZOl- If the 
polarization is normal to the free surface of the ferroelectric fllm, a geometry that 
is often used in practice, then depolarization effects must be included and this is 
the topic of the next subsection. Next we turn to epitaxial misflt strain, and how it 
can be turned to dramatically increase the spontaneous polarization of a ferroelectric 
fllm beyond that in the bulk. We end with a discussion of inhomogenous strain and 
polarization conflgurations, a more recent area of investigation that is of particular 
relevance to three-dimensional ferroelectric nanostructures |69j . 
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4.2 The Polarization at the Boundary 

The presence of the Ginzburg term, |VPp, in the Landau- Ginzburg free energy im- 
phes that a boundary condition is necessary to solve the second-order differential 
equation that results from the minimization procedure. Let us consider a semi-infinite 
ferroelectric with a second-order phase transition We begin by considering the bulk 
free energy, including the new gradient term, of this system 



■J bulk 



dV 



U I 2' ' 



4' ' 



(51) 



such that a > (implying T < Tq). We note that we can use the divergence theorem 
to split the gradient term in (j51|) into a surface and a volume integral 



dV{VPVP) = / dV V(PVP) - (PV^P) = / dS{n ■ VP)P - / dV{PV^P) 



(52) 

where n is the normal unit vector to the surface. Borrowing from analogous studies 
of superconducting magnetic [TU [72] films, we argue that there is a "surface 



tension" contribution so that the full free energy (JF = JF^f^ + J^stension) is now a sum 



of interior and surface contributions, = JF,- 
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(54) 



We note that this form of the "surface tension" free energy contribution assumes that 
there exists a temperature T*, not necessarily equal to the bulk transition temperature 
To, where the surface becomes ferroelectric; here the coefficients a and (3 are usually 
considered to be thickness- independent with (3 « a |T7t IHTj. 

If we vary the free energy density with respect to JP, we obtain 



dV 



-7V2p-a|P| +6|P|^ {6P) 



dS [7(n ■ VP) + a{T - T*)\P\\ (SP). 

(55) 

In bulk systems the surface integral in is often neglected, either due to periodic 
boundary conditions or due to the absence of variation in P on very long length-scales. 
However this approach cannot be justified for finite-size systems. 

For simplicity, let us restrict our attention here to the case where the polarization 
is only z-dependent, the surface defined as 2; = 0, and the sample occupies the space 
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(58) 



z > 0. Then the stationary condition for the total free energy to a small variation 
6P{z) leads to the second-order differential equation 

- 7—- -aP + bP^ = (56) 

together with the boundary condition 

dP 

-f— + a{T-T*)P = 0\,^o+ (57) 

at the surface. We can rewrite the boundary condition as 

dP _ P. 
dz 

which results in 

P = P(0)e"f|,^o+ (59) 

with 

6 = (60) 
a{T - T*) ^ ' 

where we note that the dimensions of 7 and a{T — T*) are and respectively 
so that [6] = L and 6 is the so-called extrapolation length; here P{0), the polarization 
at the surface, must be determined by an electrical boundary condition. We remark 
that we must obtain 

\S\ > eo (61) 

in order for this calculation to be self-consistent, where ^0 is the length-scale over 
which the polarization has spatial variation in the bulk. We see that 6 is temperature- 
dependent, and its sign is determined by the relative values of T* and T. For example 
if, as we have assumed above, T < Tq, and Tq < T*, then (T — T*) < and 6 < 0. In 
this case, the surface becomes ferroelectric at a temperature higher than that in the 
bulk, so the polarization at the surface increases as indicated in ()59|) . By contrast, 
it will decrease at the boundary if Tq > T* . This approach has been generalized 
to ferroelectric films I71H I7H] with the appropriate boundary conditions at the 
surfaces/interfaces. The solution of (jSUj) and the film analogues of leads to a 
thickness-dependent critical temperature that is qualitatively similar to that observed 
[23 HOI UBI • It can also be extended to the case where the bulk material has a first- 
order transition [Tl]; now higher-order terms in the Landau-Ginzburg free energy 
must be retained. Another generalization is to superlattices with alternating films of 
two different materials [731 IZH] • 

The scheme that we have just outlined assumes implicitly that the free energy of 
interest can be clearly separated into two distinct parts, Tinterior and T surf ace-, that 
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have the similar algebraic structures but different transition temperatures. More 
specifically, there must be a distinct separation of length-scales associated with the 
surface and the interior; this requirement translates into the inequality 

eo <\S\<L (62) 

where L >> 1 is the lateral dimension of the sample. To our knowledge, there is 
no direct experimental confirmation of the extrapolation length, though it may be 
indirectly accessible via x-ray diffraction [ZO] and far-infrared. We note that there 
is some suggestion |78^ that the Landau coefficients for thin films may have some 
thickness-dependence that is not included in this standard treatment. Another com- 
plementary approach is to assume a spatially varying transition temperature in the 
finite-size ferroelectric jZHl- This might be a productive way of incorporating known 
imperfections into the phenomenological theory. Finally we remind the reader that 
the Landau- Ginzburg approach assumes that the free energy can be expressed as 
a polynomial expansion of the average polarization; as finite-size effects become in- 
creasingly important such coarse-graining may no longer be possible. The possible 
breakdown of this free energy ansatz should be checked explicitly with first-principles 
methods for a variety of sizes and geometries. 

4.3 Depolarization Effects 

So far we have discussed finite-size planar ferroelectrics where the polarization has 
been parallel to the surfaces of interest. Now we switch to a geometry where the 
polarization is normal to the boundary, as displayed in Fig. |H| Physically we expect 
a build-up of free surface charge which, if uncompensated, results in a depolarizing 
electric field E^. It is then energetically difficult for the sample to sustain its uniform 
polarization, and phases with ferroelectric domains often result. Since, as we shall 
see shortly, the depolarization field increases with decreasing film thickness f7| ^1 
IT^ . its importance increases with reduced size and it must be incorporated in the 
phenomenological description of such finite-size ferroelectric systems. 

In practice these depolarization effects can be significantly reduced by the pres- 
ence of metal electrodes that provide charge compensation at the ferroelectric bound- 
aries. We now develop an expression for the Landau-Ginzburg free energy of a short- 
circuited ferroelectric capacitor that consists of a ferroelectric film bracketed by two 
metal electrodes (see Fig. IH)) For pedagogical simplicity, we'll consider the short- 
circuited case [Vext = 0) where the length-scale of the electrodes is significantly larger 
than that of the film [L >> d; L — > oo), which is in turn greater than the screening 
length (A) in the electrodes; succinctly our operating condition is then L >> d » X 
and a schematic of this situation is displayed in Fig. |H| We will also assume that 
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(a) 



L » d »^ 




Figure 8: a) Schematic of short-circuited ferroelectric capacitor with L >> d » X 
and accompanying plots of the spatial distribution of b) the charge Q, c) the potential 
V and d) the depolarization field E. Note that the charge distribution in b) would 
simply be ^-functions at the film-electrode interfaces for perfect electrodes; then V 
and E would simply be constant and zero. 
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the polarization gradient is negligible at the ferroelectric-electrode interfaces, namely 
that ^ = 0|2^_|_d, which implies that the inverse extrapolation length in ()58|) is zero 

(e.g. 6 — > oo). . We note that calculations in this simplified parameter regime 
yield results that are compatible with experiment fHlESI- The interested reader is 
referred to more complex treatments of the depolarization field for other parameter 
specifications [HUj . 

Let us consider the zero-potential (short-circuited) ferroelectric capacitor shown 
schematically in Fig. |Hla). The spontaneous polarization is displayed. In an idealized 
system with perfect metal plates, there would be complete charge compensation at the 
electrode-ferroelectric interface and thus no depolarization field in the film. However 
in realistic electrodes, screening charge resides within a small but nonzero distance 
from the ferroelectric boundary (cf. Fig. |Hlb). For a short-circuited ferroelectric 
capacitor, this spatial charge distribution leads to an associated voltage drop in the 
electrodes and a compensating depolarization potential across the film (cf. Fig. |Hlc). 
There is therefore an accompanying depolarization field which, in the case of uniform 
polarization, will simply scale inversely with film thickness (Fig. |HJ. Intuitively we 
expect the presence of this depolarization field to suppress ferroelectricity since it is 
antiparallel to the polarization. 

More formally we follow previous treatments jHO] in the specific limit L >> d » X 
to derive the total Landau free energy of the ferroelectric capacitor; this free energy 
has two contributions 

F = Fe + Ff (63) 

where Fe and Ff are associated with the electrodes and the film respectively. We as- 
sume that the charge density at each ferroelectric-electrode interface takes the simple 
Thomas- Fermi form p{z) = —e~^ where A is a screening length . Then, solving 
for the free energy associated with the field inside each electrode, we find that 

Fe = ^ (64) 

where eg and eo are the permittivities of the electrodes and of free space respectively. 
We note that the choice of eg, a quantity usually defined in the long- wavelength 
limit |H2j, is tricky since it enters this treatment via a boundary condition at the 
ferroelectric-electrode interface, and we refer the interested reader to two distinct 
approaches to this issue |^ IHO] where comparisons with experiment are made. 
Next we turn to the free energy in the film, which has the form 

I'd/2 r 1 ^ 

Ff = y_^^^ - -Ed{z)P{z) I dz (65) 

where Ed{z) is the depolarization field and the second term in ()65|) above represents its 
self-energy (and hence the factor of | [HZl)- Using the fact that there is no free charge 
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in the ferroelectric, we write the depolarization field in the form Ed{z) = Eq ~ 
where by continuity at the interface Eq = i?(±|) = ^. Now we need to find Eq in 
terms of P{z) since we do not know Q. Use of the short-circuit condition leads to the 
expression 

which can then be inputted into the expressions (|64j) and (|65p for Eg, and Ef respec- 
tively. We note that for P uniform and X << d, the depolarization field takes the 
form 



which varies inversely with the film thickness (d) and vanishes in the case of perfect 
electrodes (A = 0); we note that a similar result for Ed is found by introducing an air 
gap of thickness A between the film and perfect electrodes |22l IHl] thereby displacing 
the compensating charge from the ferroelectric interface. 

The Landau free energy for the ferroelectric capacitor then is 




J^pdz] + AE 
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(69) 



where we reiterate that this free energy is for the limits L » d » X where L ^ oo 
and 6 oo. We see from (I69|) that if P is uniform and A = 0, then AE = 0. However 
if P = P{z), AE is finite even for perfect electrodes [nZ|. For uniform P and finite A 
where ^ << 1, (jEHl) yields 
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(70) 



which contributes to the coefficent of the quadratic term in the overall free energy 
density. As a result, the expression for the film's transition temperature T* is 







(71) 



which leads to 



X 

eoeedao 



(72) 
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that is shifted from its bulk value (To) by a term due to the depolarization field; 
we note that when 2A ~ rf, a more complicated expression results [23 that suggests 
that there exists a minimum thickness dmin below which ferroelectricity is unstable 
fZl^llZni- For a given material, characterized by oq and Tq, this minimum thickness 
can be tuned by the choice of electrodes [201 1^ I22j whose relevant characteristics en- 
ter dmin via A and eg. However recent experiments [S3IHniEZl indicate the presence of 
ferroelectricity below this minimum cutoff, which may be consistent with theoretical 
expectations for a polydomain ferroelectric state |8Bj. Naturally the depolarization 
effects will be significantly reduced if the ferroelectric film is treated as a semicon- 
ductor jSn] rather than as an insulator. Before closing this section, we note that 
recently the dielectric properties of ferroelectric-paraelectric multilayers have also 
been studied within this type of thermodynamic approach, and it has been found 
that there exists a critical paraelectric film thickness such that ferroelectricity is com- 
pletely suppressed due to interlayer electrostatic interactions |90J. We emphasize that 
the Landau- Ginzburg approach presented to this point has assumed a single-domain 
(homogeneous) ferroelectric phase, and its generalization to include inhomogeneous 
effects will be discussed later in this chapter. 

4.4 Misfit Epitaxial Strain 

The thermodynamic behavior of a ferroelectric film is also sensitive to mechanical 
boundary conditions due to the strong coupling between the polarization and the 
strain. Application of pressure, both hydrostatic and biaxial ^Ei, is known to 
affect the transition temperature of bulk perovskite ferroelectrics. More recently ho- 
mogeneous epitaxial strain associated with substrate-film lattice mismatch has been 
included in a Landau free energy [01]; this approach results in strain-induced shifts 
in the transition temperature and the spontaneous polarization whose signs depends 
on details of specific elastic compliances and electrostrictive constants. Furthermore, 
theoretical studies [HH IHSl IHHl IHll of temperature-misfit strain phase diagrams in- 
dicate that there are observable "epitaxial phases" that would not be stable in the 
bulk. Indeed the effects of homogeneous misfit strain have been characterized and 
controlled to such a degree that BaTiO^ films have been strain-engineered so that 
their observed spontaneous polarization is more than twice that in the bulk fM^. Fur- 
thermore room-temperature ferrolelectricity has been observed in strained SrTiOs 
films j24], where we note that bulk SrTiO^ remains paraelectric down to the lowest 
observable temperatures P7] . 

Let us be more specific. In the phenomenological treatment of epitaxially strained 
ferroelectric films grown on thick substrates [HU 1021; the finite-size ferroelectric is 
approached as a bulk material with homogeneous elastic terms constrained to match 
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the substrate lattice conditions. Here it is implicit that 



L» dc> d> (73) 

where L and d are the substrate and the film thicknesses respectively, is the corre- 
lation length and dc is the critical film thickness above which elastic defects appear; 
dc varies roughly as the inverse of the lattice mismatch and is usually approximately 
lOnm pn]. With this set of modelling assumptions, the in-plane strains, r^i, ri2 and 
rjQ (here we follow the conventional Voigt notation) are constant throughout the film 
thickness and are completely controlled by the substrate-film lattice mismatch. The 
associated stresses ui, (T2 and are finite but not fixed. By contrast, 

= ^4 = as = (74) 

since there are no tractions acting on the top film surface. For the special case of a 
(001) ferroelectric film grown on a cubic paraelectric substrate 

a — ao 

111 = 112 = 11 = (75) 

a 

where a and ao are the constrained and the free film lattice constants respectively. 
For this simple case of cubic symmetry 7/6 = since the angle between the two 
lattice vectors remains unchanged {0 = f ). (f7^ and (ffH|l represent mixed mechanical 
boundary conditions associated with two-dimensional clamping, and the standard 
elastic free energy J-'{P, cr) cannot be used to find the equilibrium properties of these 
systems I91j. Instead a Legendre transformation to a modified thermodynamic 
potential 

T = T \ r]iai + 7720-2 + r/ecre (76) 

must be performed in order to study the equilibrium properties of this constrained 
film. 

For pedagogical simplicity, we consider a uniaxial ferroelectric where P is the 
polarization in the z-direction. The free energy, with condition (f7^ . of a cubic ferro- 
electric is 191 



= J^P- ^Sii(o-i + al) - Ql2{(cri + a2)P^} - Sl2Cri0'2 - ^S440-g (77) 



where Qij and Sij are the electrostrictive constants and the elastic compliances at 

dcTi 



constant polarization respectively. Using jf- = —T]i, and solving for ai = a2 = cr (in 



this special case cxe = 0), we find that 



^ = + l~aP' + hp' + IcP' (78) 

S11 + S12 2 4 6 ^ ^ 
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where 

a = a (79) 

Sll + Si2 

and 

b = b+ ^^'^ (80) 

Sll + Su 

so that the coefficients of both the quadratic and the quartic polarization terms 
in the Landau free energy are renormahzed. This has important imphcations for 
the thermodynamic properties of the thin ffim; more specifically it means that its 
transition temperature is shifted from the bulk value (Tq) 

T* = To+ /"^^'^ , (81) 
a{Su + S12) 

as is its spontaneous polarization (please see the discussion preceding (fTUIl) due to 
the in-plane straining of the ffim by the compressive substrate. The signs of these 
shifts will be determined by the relative signs and magnitudes of Qu, sn and S12. 
We note that such misffi epitaxial strain can change the nature of the transition from 
first- to second-order (6 < but b > 0) For ferroelectric films with multicomponent 
polarization and with different orientations, this approach yields a temperature-strain 
phase diagram with equilibrium phases that are not stable in a mechanically-free bulk 
sample [HIl IHS UHl IMl US] • The resulting low-temperature phase diagrams are very 
sensitive to the Landau coefficients, particularly for increasing misffi strain; here the 
phenomenological approach is nicely complemented by ffist-principles studies that 
can resolve uncertainties associated with these input parameters [94j. This combined 
phenomenological-first principles approach has been very successful in explaining the 
phase behavior of epitaxially strained films and we refer the interested reader to 
a sample of recent experimental papers in this field jHll 1111 I100| IIOH llU2j . For 
ultrathin ferroelectric ffims with biaxial compressive strain at the substrate interface, 
the competition between elastic and surface effects has also been considered ||103] . In 
particular we note that if we do include a gradient (Ginzburg) term to the free energy 
then the epitaxial strain leads to a modified correlation length 

(82) 

where ^ = is that associated with the bulk material that is mechanically free |lU3j . 
Similar mechanical boundary conditions have also been studied in epitaxially-strained 
superlattices |lU4j . We note that recently a combined numerical-thermodynamic ap- 
proach resulted in generalized temperature-strain phase diagrams for thin ferroelectric 
ffims that include multidomain phases |lU5j inaccessible by the approaches described 
here that assume homogenous polarization. 



35 



4.5 Inhomogeneous Effects 

In the previous discussion, we assume that the energy cost of uniformly straining 
the film to lattice-match the substrate is significantly less than that associated with 
other strain relaxation mechanisms; these include the creation of elastic dislocations, 
polydomain formation and even multiphase coexistence which are all inhomogeneous 
in nature. We expect elastic homogeneity for films thinner than a critical thickess 
{d < dc), where d^ scales (roughly) inversely with the lattice mismatch Typ- 
ically such coherent epitaxially strained films have thicknesses c/ ~ 10 nm, though 
film growth parameters can be tuned to achieve ci ~ 50 nm by kinetically suppressing 
other elastic relaxation processes inHl- However the physics of thicker ferroelectric 
films, particularly those used in current applications {d ~ 120 nm) [58j, does include 
inhomogenous effects and next we review aspects of this topic within a phenomeno- 
logical framework. 

The misfit epitaxial strain decreases roughly exponentially |1()6| ll()7[ ll()8j with 
distance (fjiz) ~ fjoe'^) from the film-substrate interface with a "strain depth" of 
approximately A ~ 300 nm and thus is negligible for d » X. In films of thickness 
dc < d < X, elastic defects often form to accomodate film-substrate lattice mismatch. 
While permitting the film's lattice constants to relax, these defects do generate in- 
homogenous strains which couple to the film's polarization and therefore affect its 
ferroelectric properties. We note that compositional and thermal gradients can also 
produce such inhomogeneities |109j . Here we will focus on such elastic effects. Indeed 
the epitaxial phenomenological treatment can be generalized to include isotropic inho- 
mogenous strain due to such lattice defects |ll()j : this results in further contributions 
to the coefficients of both and in the free energy, and thus to a defect-induced 
change in the transition temperature qualitatively consistent with experiment |110j . 

The broadening of the dielectric peak in ferroelectric films with decreasing thick- 
ness has been observed by several groups |108t IllH I112| I113j . In bulk systems, 
such behavior is often attributed to disorder and/or low-dimensional effects jI14] but 
these explanations are not appropriate to ferroelectric films due to the underlying 
long-range interactions. Indeed we expect these dipolar systems to be more sensitive 
to electromechanical boundary conditions than to local fluctuations. This qualitative 
idea has been confirmed by measurements on ferroelectric free-standing lamellae of 
sharp dielectric peaks, strongly indicating that it is due to substrate-related effects 
most probably associated with interfacial lattice mismatch |115j . 

Recent experiments on the fiexoelectric coupling between strain gradients and the 
polarization indicate that these effects could be important in thin films |116t IllTj . 
Theoretically it has been shown that this fiexoelectric effect is enhanced in high- 
permittivity materials such as ferroelectrics jll8j where the coupling is maximized 
near the dielectric peak. This fiexoelectric coupling has been incorporated into a 
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Landau- Ginzburg framework with the goal of studying its effect on the thermody- 
namic properties of epitaxially strained ferroelectric films |119j . The strain gradient 
contribution, /S.fc = 7^^^; couples to the polarization in the free energy, acting as 
an effective field, so that the equation for the spontaneous polarization 

~aP + hP'' = Afc\p=Po (83) 

can only have finite solutions, even in the absence of an external electric field (i? = 0); 
here a and h are the Landau coefficients that have been renormalized by epitaxial 
strain. As a result the inverse dielectric susceptibility 

X"' = 36Po' + « (84) 

can never be zero, so the standard singularity in x is replaced by a broadened peak. 
Similarly the temperature scale associated with the onset of reversible polarization is 
distinct from that of the maximum of x- This approach also yields a polarization that 
increases with decreasing film thickness despite the associated reduction and broaden- 
ing of the dielectric peak; these seemingly contradictory features are consistent with 
known experiment |119j . 

We have just seen that the presence of strain gradients via the flexoelectric effect 
leads to the broadening of the dielectric peak in ferroelectric films. This phenomena 
can be generalized to consider other effective fields that couple linearly to the po- 
larization |12Uj : possible origins of this effective field include asymmetric electrodes, 
compositional and temperature gradients as well as the stress profiles already con- 
sidered. Qualitatively the effects on the dielectric response will be the same, though 
naturally quantitative differences will depend on the physical origin of the effective 
field. 

Inhomogeneities, particularly strain and temperature gradients, can also lead to 
multidomain formation |122j as an overall energy reduction mechanism; this behavior 
has been observed in relatively thick {d ~ 500 nm) films |123j . The competition 
between different strain relaxation processes including domain formation and mis- 
fit dislocatio has been studied theoretically |124[ I125j . Within the framework of a 
Landau- Ginzburg theory, the development of domain structure has been addressed 
|88| I12H I122j ; qualitatively this is done by seeking polarization solutions of the form 

P = P^ + 5Pk{z)e'^^ (85) 

where 5Pk{z) « Pq. Such a treatment in nearly cubic ferroelectrics indicates that a 
multidomain state may well be stable down to atomic film thicknesses, length-scales 
well below the previous estimates that only described single-domain ferroelectric- 
ity [HH]-The time evolution of the polarization and therefore the resulting domain 
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structures can also be obtained by numerical studies of the time- dependent Landau- 
Ginzburg equations[125]; such phase-field studies predict specific domain morpholo- 
gies jSni EEl I127| Il28j for epitaxial thin films of a variety of materials including 
PbZri-xTixOs, SrTiO^ and BaTiO^. We note that inhomogeneities in multicom- 
ponent ferroelectrics may also lead to multiphase coexistence [125. .129^ .130. .131j . 
particularly in the region of a first-order phase transition. 

5 Summary and (Some) Open Questions 

At this point we have discussed the main features associated with ferroelectrics that 
have been addressed within a phenomenological Landau- Ginzburg theory. Due to 
their underlying long-range interactions, we have learned that bulk ferroelectrics are 
well described within this framework; for the same reason, their finite-size counter- 
parts are very sensitive to electromechanical boundary conditions including surface 
charge and misfit epitaxial strain. Inhomogeneous effects, particularly strain gradi- 
ents, have also been discussed. We have now worked our way to the present and 
naturally there remain a number of challenges that can be addressed within this 
framework. Let us discuss a few of these: 

(i) Strain-gradients and finite-size effects in novel nanogeometries 

Three-dimensional ferroelectric nanogeometries are being explored as key compo- 
nents in competitive high-density data storage devices |132j . The complexity of their 
topologies and their boundary conditions tends to favor inhomogenous polarization 
and strain configurations (e.g. |133| I134j ) which remain largely unstudied, partic- 
ularly from a phenomenological standpoint. For example, the energetics associated 
with the formation of misfit dislocations in ferroelectrics with curved topologies is 
known to be different than those in planar structures |135j , and the associated strain 
gradients could lead to complex three-dimensional polarization patterns. 

(ii) Multiferroics 

There exist a number of materials that display both long-range ferroelectricity 
and incommensurate magnetic order |136j and we refer the interested reader to the 
chapter by N. Spaldin on this topic. Landau treatments |137| I138j for this class 
of materials have provided symmetry-based arguments for the type of underlying 
interactions that must be present at the microscopic level; furthermore the couplings 
between ferroelectric and magnetic order parameters can be studied in nanostructured 
electromagnets |139j . Detailed interplay between microscopic and phenomenological 
models should provide guidance on material parameters necessary to enhance the 
magnetoelectric couplings, particularly at room temperature for device applications. 

(iii) Dynamics 

The dynamical dielectric response of ferroelectrics should be accessible via time- 
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dependent Landau-Ginzburg theory, similar to studies performed with this method in 
ferroelastics (e.g. [3T]). More generally a detailed phenomenological study of domain 
nucleation |14nj . domain motion (and domain interlocking J141]) could be particularly 
useful for the characterization of switching properties of thin-films [ 142j . Also glassy 
behavior has been observed in a number of ferroelectric relaxor materials |143j . and 
a phenomenological model of these relaxational dynamics should be feasible. 

(iv) Possible Breakdown of Landau-Ginzburg theory 
Landau-Ginzburg theory may well break down in ultrasmall ferroelectric struc- 
tures where the necessary averaging of the order parameter (i.e. the polarization) is 
no longer possible. Also it remains a challenge to see how this highly local approach 
can be used to described nonlocal elastic and Coulomb effects, particularly vis a vis 
domain formation and energetics. 

(v) Quantum Critical Fluctuations 

The paralelectric-ferroelectric transition in the limit T — may well be a textbook 
candidate |144U145|ll4tij for quantum critical behavior |147j . Though some theoretical 
predictions |148j exist, logarithmic corrections to mean-field (Landau) scaling should 
be both theoretically and experimentally accessible and have yet to be explored. As 
in the classical case ferroelectric materials may well provide the setting for the 
first detailed study of criticality at a quantum phase transition. 

These are only some of the many outstanding challenges that the Landau-Ginzburg 
theory of ferroelectrics still has to address. Most generally, it can provide a solid bridge 
between the physics on atomistic scales and macroscopic measurable quantities, par- 
ticularly in finite-size ferroelectrics where electromechanical boundary conditions are 
crucial and therefore the physics spans many length-scales. In closing, we invite the 
(hopefully still) curious reader to partake in these and other adventures with this 
phenomenological approach to ferroelectric materials. 
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